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Abstract 

We propose a new method to identify transitions from a topological insulator to a band insulator 
in silicene (the silicon equivalent of graphene) in the presence of perpendicular magnetic and 
electric fields, by using the Renyi-Wehrl entropy of the quantum state in phase space. Electron- 
hole entropies display an inversion/crossing behavior at the charge neutrality point for any Landau 
level, and the combined entropy of particles plus holes turns out to be maximum at this critical 
point. The result is interpreted in terms of delocalization of the quantum state in phase space. 
The entropic description presented in this work will be valid in general 2D gapped Dirac materials, 
with a strong intrinsic spin-orbit interaction, isoestructural with silicene. 
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I. INTRODUCTION 


Silicene is a two dimensional crystal of silicon which has been studied theoretically OEI 
and recently experimentally [3H7]. The low energy electronic properties can be described by 
a Dirac Hamiltonian, like in graphene, but the electrons are massive due to a relative large 
spin-orbit coupling Ago- In fact, silicene has a band gap | Ag^| (s and ^ denote spin and valley, 
respectively) which can be controlled by applying a perpendicular electric held £z = 

(/ is the interlattice distance of the buckled honeycomb structure) to the silicene sheet (see 
[^for a plot of the band gap). It has been demonstrated |H] that there is a quantum phase 
transition, from a topological insulator (TI, | A^l < Ago) to a band insulator (BI, |A^| > Ago), 
at a charge neutrality point (CNP) aI^^ = s^Ago, where there is a gap cancellation between 
the perpendicular electric held and the spin-orbit coupling, thus exhibiting a semimetal 
behavior. 

A 2D topological insulator is also known as a quantum spin Hall state, and it was the¬ 
oretically studied in [H] and hrst discovered in HgTe quantum wells in uni. The common 
characteristic of a TI-BI transition is a band inversion with a level crossing at some critical 
value of a control parameter (electric held, quantum well thickness, etc). In this paper we 
hnd that electron-hole entropies of a quantum state (in a phase-space representation) also 
exhibit this crossing/inversion behavior at the CNP, thus providing a new characterization 
of the TI-BI quantum phase transition (QPT). 

The study of phase space properties in quantum systems can be done by using the Husimi 
function of a quantum state \'ip), which is dehned as the squared overlap between a minimal 
uncertainty (coherent) state |a) and IV’), thus giving the probability of hnding the quantum 
system in a coherent state. Husimi quasiprobability distribution has been useful for a phase- 
space analysis of metal-insulator transition HU, for the study of quantum chaos in physics 
j, and recently in the analysis of QPTs in algebraic models [TBtiTB] . Renyi-Wehrl entropies 


of the Husimi function (see section HI for a dehnition) have been proved to be a good marker 
for QPTs of some paradigmatic models like: Dicke model [13] of atom-held interactions 
(undergoing a QPT from normal to superradiant), vibron model [H] for triatomic molecules 
(undergoing a shape QPT from linear to bent) and Lipkin-Meshkov-Glick model [T5| used 
for example in nuclear physics and quantum optics. Renyi-Wehrl entropy displays either 
an extreme, or an abrupt change, at the critical point and therefore provides a detector of 
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QPTs across the phase diagram of a critical system. 

In this work we shall analyze the topological-band insulator transition in silicene by 
using these phase-space techniques. The paper is organized as follows. Firstly, in Section 
[m we shall remind the low energy Hamiltonian for silicene, its eigenvectors and eigenvalues. 


Then, in Section III, we will dehne the Husimi function in this model and the Renyi-Wehrl 


entropies of the Husimi function. In Section W we will show that the Renyi-Wehrl entropy 
turns out to be a clear marker of topological-band insulator transitions in silicene. Finally, 
some concluding remarks will be given in the last Section. 


II. LOW ENERGY HAMILTONIAN 


Let us consider a monolayer silicene with external magnetic B and electric 8^ helds applied 
perpendicular to the sheet. The effective low energy Hamiltonian is given by [S] 

Hi = VF{(yxPx - if^yPy) “ ^sAgoCr^ + (1) 


where ^ corresponds to the inequivalent corners K {^ = 1) and K' (^ = — 1) of the Brillouin 
zone, respectively, aj are the usual Pauli matrices, vp = 5 x 10^ m/s is the Fermi velocity of 
the Dirac fermions (with up and down spin values being represented by s = ±1, respectively) 
and Aso is the band gap induced by intrinsic spin-orbit interaction, which provides a mass to 
the Dirac fermions. The spin-orbit energy gap induced by intrinsic spin-orbit coupling has 
been estimated (using density functional theory and tight-binding calculations piTHTU] ! as 
Ago ~ 1.55 — 7.9 meV, and we will consider Ago = 4 meV along this paper. The application of 
an external constant electric held 8z creates an potential difference Ax = ISz, with I ~ 0.23A, 
between sublattices yielding a tunable band gap. See Figure [^for a graphical representation. 

Using the Landau gauge, A = (0, Bx, 0), and doing the minimal coupling p —)■ p -f for 

the momentum, one can easily diagonalize the Hamiltonian ([^ with eigenvalues [H |20] 


= 


sgn(n) 




n ^ 0, 


( 2 ) 


[ -^Ag^, n = 0, 

in terms of the Landau level index n = 0,±1,±2,..., the cyclotron frequency u = 
VF\/2eB/h and the lowest band gap A^^ = (A^ —s^Ago)/2. The corresponding eigenvectors 
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FIG. 1: The band gap = (A^ — s^Aso)/2 as a function of the electric potential A^ created 
by the perpendicular electric field applied to the silicene sheet. Ago is the spin-orbit energy gap 
induced by the intrinsic spin-orbit. 


for the and K' points are 


f-afllnl-(1+0/2) 
B;«||n|-(1-0/2) 


(3) 


where the constants and are given by 

245 = 


n^O, 


s5 = 


(1-0/2, n = 0, 

|E°^|-sgn(n)Aag 


2|Kh 

(1 + 0/2, n = 0. 


, n 7^ 0, 


(4) 


and ||?7,|) denotes an orthonormal Fock state of the harmonic oscillator. 

In Fignrej^we represent the low energy spectra given by eq. ([^, as a fnnction of the 
external electric potential for B = 0.5 T (where T is the nnit symbol for Tesla). The 
transitions TI-BI occnr at the critical points = s^Ago, where there is a band inversion 
for the n = 0 Landan level (either for spin np and down) at both valleys. The energies 
and Eq^’^ have the same sign in the BI phase (|A^| > Ago, in bine), and different sign in the 
TI phase ( |A^| < Ago, in pink). 


III. INVERSE PARTICIPATION RATIO AND RENYI-WEHRL ENTROPY 


The standard canonical (or Glanber) coherent states for the ID qnantnm harmonic os¬ 
cillator are given by 

OO n 

|a) =e-l“lV2gaat|Q^ « 1^^^ (5) 

n=0 
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FIG. 2: Low energy spectra of silicene as a function of the external electric potential for B = 0.5 
T. Landau levels n = ±1,±2 and ±3, at valley ^ = 1, are represented by dashed (electrons) and 
solid (holes) thin lines, black for s = — 1 and red for s = 1 (for the other valley we simply have 
En~^ = The lowest Landau level n = 0 is represented by thick lines at both valleys: solid 

at .^ = 1 and dashed at ^ = — 1. Vertical blue dotted grid lines indicate the CNPs separating BI 
(blue) from TI (pink) phases. 

with a G C (the phase space). The coherent states form an overcomplete set of the corre¬ 
sponding Hilbert space and verify the closure relation: 

1= [ |a)(a|—, (6) 

with (Pa = (iRe(Q:)dIm(a). Using the coherent state (Bargmann) representation, we can 
associate a quasiprobability distribution (the so called Husimi function) (5p(a) = to 

every normalized harmonic oscillator state ip). This dehnition is straightforwardly extended 
to any density matrix p as Qp{a) = (a|p|Q:). For example, the Husimi function of a Fock 
state \n) is simply 

<J„(a) = |{n|a)P=U)!-|ap“ (7) 

n\ 

Coherent states are said to be “quasi-classical” because of their minimum uncertainty (and 
dynamical) properties. To quantify the spread/localization of a wave packet p in phase 
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( 8 ) 


space, the z/-th moments of the corresponding Husimi function Qp, dehned by 

m; = f —(<^♦(a))^ 

are often used. Note that, using (|^, we have = 1 for normalized ip. We shall focus on 
the u = 2 moment, also called the inverse participation ratio, which is related to the inverse 
area in phase space occupied by the Husimi function. Namely, for a Fock state \n) we have 


Mp 


(2n)! 

22n+l(^!)2’ 


(9) 


which is maximum for n = 0, Mg = 1/2 (minimum area). Instead of the momentum M^, 
we shall use the Renyi-Wehrl entropy W^, = for the sake of convenience. Note 

that, for 1 / —)■ 1, the previous expression reads = —f Qjf,(a) In (Q^(a)) d^a/7r, the so 
called Wehrl entropy. As conjectured by Wehrl [2l] and proved by Lieb [22], any Glauber 
coherent state la') has a minimum Wehrl entropy of IF/, = 1. We shall use u = 2, for which 
the minimum Renyi-Wehrl entropy is W/iin = ln(2), corresponding to a coherent state. 


IV. RENYI-WEHRL ENTROPY AS A MARKER OF TOPOLOGICAL-BAND IN¬ 
SULATOR TRANSITIONS 


Now, we shall calculate the Renyi-Wehrl entropy IF/ of a Hamiltonian eigenstate (|^ as 
a function of the electric potential in order to analyze the TI-BI transition. We shall 
restrict ourselves to the valley ,^ = 15 omitting this index from ([^ and (|^. All the results 
for ^ = 1 are straightforwardly translated to ^ = — 1 by swapping electrons for holes (i.e., 
n 0 —n) and spin up for down (i.e., s <H- —s). Taking into account equations ([^ and ([^, 
we can write the expression of the Husimi function of a Hamiltonian eigenstate lu)^ as 

Qp{a) = \{a\n)s\‘^ = (A/)^QH-i(a) + (R/)2Q|n|(a), (10) 

and the second moment as 


M/’" = 


(A/)^MJ|_, + (i?:)2M|/| 


+2{AX 


(2( 

n\ 

— 

1)! 

4kl 

\n\ 

|!(|n| 

-1)! 


( 11 ) 


where we are using ([^ and the value of J ^Q\n\{oi)Q\n\-i{oi)■ For the sake of convenience, 
we shall use the Renyi-Wehrl entropy IF/’'^ = — ln(M/’*) to make graphical representations 
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FIG. 3: Renyi-Wehrl entropy Wn^ of a Hamiltonian eigenstate |n)sj as a function of the electric 
potential for the Landau levels: n = 1, 2 and 3 (electrons; dashed curves) and n = — 1, —2 and 
—3 (holes; solid curves) and valley ^ = 1. Black and red curves correspond to spin s = — 1 and 
s = 1, respectively. Electron and hole entropy curves cross (“entropy inversion”) at the critical 
value of the electric potential = —sAgo (vertical blue dotted grid lines indicate this CNPs), 
with crossing entropy values denoted by Cn = W’n’*(Ai°^), n = 1,2,3 and given in equation © 
(horizontal blue dotted grid lines indicate these crossing points). Thin and thick lines correspond 
to magnetic fields B = 10“^ and B = 10“^ T, respectively. 


as a function of for different values of the Landau level n, spin s = ±1 and the magnetic 
held B. In Figure we plot as a function of A 2 for n = 1,2,3 (electrons) and 

n = — 1, —2, —3 (holes) with spin up s = 1 (in red) and down s = — 1 (in black). Entropy 
for holes (n < 0, solid curves) is an increasing function of the electric potential A^, whereas 
entropy for electrons (n > 0, dashed curves) is a decreasing function of A^. Electron-hole 
entropy curves cross at the CNP Ai°^ = —sAgo, where the entropy is 

’(8|n| — 3)r(|n| — I 




16^/nT{\n\ -b 1) 


( 12 ) 
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FIG. 4: a) Same representation as in Figure]^ but for n = ±1 only and magnetic fields B = 0.1 
(thin) and B = 1 T (thick), respectively, b) Combined entropy Wn^ + of electrons plus holes 
for n = 1. The entropy of the electron-hole compound exhibits a maximum of 2Cn at the CNF 
— -sA 

LAz — *^SO* 


Note that the critical entropy (12) only depends on the Landau level n (and not on other 
physical magnitudes like: B,VF,e, s, etc) thus providing a universal characterization of the 
topological-band insulator QPT, which is shared with other two-dimensional crystals like 
Ge, Sn and Pb. Indeed, in Figure|^we plot entropy for magnetic held B = 10“^ (thick lines) 
and B = 10“^ T (thin lines) and note that entropy slopes at decrease with B but they 


have a common critical entropy (12). For electrons (holes), minimum (maximum) entropies 


for a given Landau level n are attained at asymptotic values —)■ cx), and maximum 

(minimum) entropies are attained at A^ —)• —oo. The minimum and maximum entropies 





















turn out to depend only on \n\ as 


W, 


2,s ^ f l^|r(|n| - 1) 

V2y^r(|n| + i) 


ty2,s = _ In 
n,max 


(|n| - |)r(|n| - I) 

2^/TTT{\n\ + 1) 


( 13 ) 


t 2 ,s 


n+l,min = aS Can alsO 


where T denotes the usual gamma function. Note that 
be seen in Figure js For example, = ln(2) (which means that the Hamiltonian 

eigenstate | ± 1)^ is a coherent state of electrons for —)■ cx) and a coherent state of holes 

for — >■ —cx)) and hF^fmax = 21n(2) (which means that | ± 1)* is a “cat state” of electrons 
for A 2 —)■ —cx and a cat state of holes for A^ —)> x). Indeed, as commented after (|^, 
the minimum entropy of = ln(2) is attained for coherent states. The so-called “cat 

states” arise in many models undergoing a QPT (see e.g., [I3HIS]) and are a superposition 
of two coherent (semiclassical) states with negligible overlap, so that their entropy is twice 
the minimum entropy, that is, 21n(2). 

In Figure]^ we plot the case n = ±1 for higher magnetic fields (H = 0.1 and H = 1 T in 
panel a) and the combined entropy of electrons plus holes, IF^’^ -f IF"^’,^ (in panel b), which 
exhibits a maximum of 2C'„ at the CNPs. For n = ±1, the maximum entropy of the electron- 
hole compound is 2Ci = 21n(16/5) ~ 2.3263 (see Figure]^ panel b). The asymptotic 
behavior of the maximum combined entropy for large n is 2C'„ ~ ln(47r)-|-ln(|n|)-|-0(l/|?7,|)^/^, 
which shows an increasing logarithmic behavior with |n|. 

The existence of a maximum electron plus hole entropy means that the combined wave 
function spreads across the phase space, taking up a bigger volume (it is more delocalized) 
at the CNP. In Figure we also plot the combined renormalized entropy 


wZ = + if!’: - - IF, 


7-2,5 

n,min 


2,5 

n,max 


(14) 


of particles plus holes for n = 1 (solid lines) and n = 2 (dashed lines) together. We see that 
entropy “hats” become flatter and flatter as n increases. 

The global behavior of the entropies IF^’^, as a function of the electric held strength, 
exhibits some common features with the band inversion behavior for the n = 0 Landau 
level (see Fig. |^and comments at the end of Section 0, Indeed, dehning the renormalized 
entropy 


IF:’^(A,) = IF:’^(A,)-a = -ln 


M!’-(A , 

M!’^(Ai 


(15) 
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FIG. 5: Renormalized Renyi-Wehrl entropies (14) of the electron-hole compound for a Hamiltonian 
eigenstate as a function of the electric potential (R = 1 T) for the Landau levels n = 1 

(solid lines) and n = 2 (dashed lines), and valley ^ = 1. Black and red curves correspond to spin 
s = — 1 and s = 1, respectively. Vertical blue lines indicate the CNPs. 


one can see that the quantities and (for n = ±1,±2,...) have the 

same sign in the BI phase (jA^I > Ago), and different sign in the TI phase (|A^| < Ago). See 
e.g. Figure]^ for i? = 1 T (thick lines). Therefore, this “entropy-inversion” provides a clear 
characterization of each phase. 

Additionally, taking into account that measures the area occupied by the eigenstate 
|n)g in phase space, we can say that, in the BI phase, an spreading (resp. localization) 
effect can be observed in the electrons (resp. holes) for A^ < —Ago and in the holes (resp. 
electrons) for A^ > Ago; in the TI phase there is a spreading effect in the electrons with spin 
up and holes with spin down, and a localization in electrons with spin down and holes with 
spin up. We observe that the lower the magnetic held strength, the sharper these spreading 
and localization effects are. 


V. CONCLUSIONS 

In this work we have shown that the Husimi function is an appropriate tool to study 
topological-band insulator transitions in silicene and other gapped isoestrucural Dirac ma¬ 
terials with intrinsic spin-orbit coupling as germanene, stantene or Pb. In particular, the 
Renyi-Wehrl entropy of the Husimi function of order two (related to the inverse participa¬ 
tion ratio) has been used to analyze the TI to BI quantum phase transition, showing the 
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critical points of this transition clearly as a crossing/inversion behavior of the electron and 


hole entropies at the charge neutrality point. Critical (12) and asymptotic (13) values of 
electron-hole entropies only depend on the Landau level n, and not on any other physical 
quantities like: magnetic strength, Fermi velocity, electron charge, etc. Therefore, these en- 
tropic measures provide a universal characterization of the topological-band insulator QPT, 
which is shared with other two-dimensional crystals like Ge, Sn and Pb. 

We want to emphasize that all the descriptions in terms of Renyi-Wehrl entropies that we 
have presented in this work will be valid in (isoestructural) 2D gapped Dirac materials with 
a strong intrinsic spin-orbit interaction. Actually, all the characterizations studied here will 
apply to Ge, Sn and Pb counterparts. In fact, all the Figures will be valid by simply replacing 
the value of Ago in each case: = 11.8 meV for germanene, A|^ = 36.0 meV for stanene, 

A^J’ = 207.3 meV for Pb [23]. Moreover, Dirac equation (and modihcations) provides 
an effective model for a large family of topological insulators, and an entropic description 
like the one done in this letter could also sheed new light on the better understanding of 
topological phases in these physical systems. This deserves a separated and more thorough 
study and it will be the subject of future work. 
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